Critical temperature of a Rashba spin-orbit coupled Bose gas in harmonic traps 
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We investigate theoretically Bose-Einstein condensation of an ideal, trapped Bose gas in the 
presence of Rashba spin-orbit coupling. Analytic results for the critical temperature and condensate 
fraction are derived, based on a semi-classical approximation to the single-particle energy spectrum 
and density of states, and are compared with exact results obtained by explicitly summing discrete 
energy levels for small number of particles. We find a significant decrease of the critical temperature 
and of the condensate fraction due to a finite spin-orbit coupling. For large coupling strength and 
finite number of particles A*', the critical temperature scales as N'^^^ and A^^^ in three and two 
dimensions, respectively, contrasted to the predictions of A^^^'' and N^^"^ in the absence of spin-orbit 
coupling. Finite size corrections in three dimensions are also discussed. 
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I. INTRODUCTION 

The recent experiment on spin-orbit (SO) coupled 
spinor Bose gases of ^^Rb atoms |l| has stimulated great 
interest in the theoretical study of SO physics in both 
Bose-Einstein condensation (BEG) and fermionic super- 
fluidity. It is well-known that the SO coupling leads 
to many interesting phenomena in condensed matter 
physics. A typical example is the recently discovered 
topological insulators or quantum spin Hall states 0, [3 . 
In degenerate atomic gases, due to unprecedented con- 
trollability in the interatomic interaction, geometry and 
purity 0)01) the SO coupling may bring even more in- 
triguing states of matter (6l-[24|. 

For a SO coupled BEG, non-trivial structures, such 
as the jiensity-stripe state |9|. | Tl| . half-quantum vortex 
state [10] and lattice state [la. Il3|. are predicted. For 
an atomic Fermi gas near Feshbach resonances, new two- 
fermion bound states with anisotropic mass are formed 
even at a negative s-wave scattering length [l^ [2l], [S^l , 
leading to the prospect of anisotropic superfluidity with 
mixed s- and p-wave components [22|. By imposing an 
external Zeeman field, novel topological superfluid sup- 

P'ting zero-energy Majorana modes may also emerge 
) [13 nil- To observe these new states of matter, it 
is necessarily to cool the temperature below a threshold, 
which may depend critically on the SO coupling. The 
purpose of this work is to determine the critical temper- 
ature of trapped atomic Bose gases with Rashba type 
SO coupling. We focus on an ideal, non-interacting Bose 
gas, since the critical temperature is less affected by weak 
interatomic interactions |25| . 

Theoretically, the critical temperature of a homoge- 
neous Bose gas is greatly suppressed by the Rashba 
SO coupling, as the low-energy density of states (DOS) 
is dramatically modified jH, ll8| . In three dimensions 
(3D) without Rashba SO coupling, the low-energy DOS 
p{E) vanishes as \/E. As a result, the number of to- 
tal particles occupied at finite energy levels, given by 
N{T) = dEp{E)/{e^/''^^ - 1), saturates at finite 



temperature T j^]. This leads to the well-known macro- 
scopic occupation of the ground state, i.e., the formation 
of a BEG. In the presence of Rashba SO coupling, how- 
ever, the low-energy DOS becomes a constant (see Ap- 
pendix A) [^, [3 ) reminiscent of a two-dimensional (2D) 
system. The thermal occupation N(T) can be logarith- 
mically divergent. The critical temperature is therefore 
precisely zero, ruling out the possibility of BEG at any 
finite temperatures [4j. 

In this paper, we show that in the presence of a har- 
monic trap the Rashba SO coupling does not destroy the 
BEG at finite temperatures, as the thermal occupation 
N{T) remains finite. Actually, the critical temperature 
is not affected by the Rashba SO coupling in the ther- 
modynamic limit where the number of particles TV be- 
comes infinitely large. This is because the occupation of 
low-energy states, modified by the SO coupling, becomes 
negligible as — > oo. However, in the experimentally 
relevant situation in which numbers of particles range 
from a few thousand to a few million, we find a sig- 
nificant decrease of the critical temperature and of the 
condensate fraction. In particular, at a sufficiently large 
Rashba SO coupling strength, the critical temperature 
scales like N^/^ and N^^^ in three and two dimensions, 
respectively, in sharp contrast to the scaling of N^^^ and 
A^i/2 without SO couphng ^i, |26h28|]. We derive these 
results either by summing discrete energy levels for small 
number of particles or by using a continuous DOS un- 
der the semi-classical assumption that the level spacing 
is negligible compared to the temperature. The former 
approach also enables the investigation of finite-size cor- 
rection to the critical temperature. 

The paper is structured as follows. In the next section 
(Sec. II), we introduce the theoretical model for Rashba 
spin-orbit coupled ideal Bose gases in harmonic traps and 
solve the single-particle energy spectrum. In Sec. HI, we 
present the 2D and 3D DOS with or without the contin- 
uous spectrum approximation. The critical temperature 
and condensate fraction are then calculated in Sec. IV 
for both 2D and 3D cases. Next, the finite size effect in 
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3D is discussed in Sec. V. Finally, Sec. VI is devoted 
to conclusions. The calculation of the DOS of a homo- 
geneous 3D Rashba SO coupled system is given in the 
Appendix A. 



II. MODEL HAMILTONIAN AND 
SINGLE-PARTICLE ENERGY SPECTRUM 

We consider a two-component (spin-1/2) Bose gas in 
2D or 3D harmonic traps, V2D{r±) = Mujj_{x'^ +y'^)/2 = 
Mio\r\/2 or Vzoir^.z) = M(cj^r^ +^^2:2^2, with a 
Rashba SO coupling Vso — —i^Ris'xdy — ^ydx) in the 
x~y plane, where A/j is the Rashba SO coupling strength 
and &xi o'y, and az are the 2x2 Pauli matrices for pseudo- 
spin. The model Hamiltonian for single-particle is de- 
scribed by. 



ns = 



-h^V'^l2M + Vt -iXnidy + id^) 
-iXuidy ~ idx) ^tf\/^/2M + Vt 



(1) 



where the trapping potential VT(rj_) = V2D{r±) in 2D 
and VriT) — V3]j{r±, z) in 3D harmonic traps. The char- 
acteristic length scales of harmonic traps in the x — y 
plane and z-direction are given by, a±_ = \J h/ {Mu}±_) 
and az = \Jh/{Mujz)^ respectively. For the SO cou- 
pling, we take a dimensionless coupling strength A50 = 

In the 2D case, it is convenient to use polar coordinates 
= (rj_,(/7). in which —i{dy ± idx) — e^'^'^[±d/dr± — 
(i/rj_)d/d(p]. As the harmonic potential is isotropic, the 
single-particle wave- function has a well-defined azimuthal 
angular momentum 1^ = m and takes the form. 



0m (r±) 



'2-k' 



(2) 



which preserves the total angular momentum jz = Iz + 
Sz = m + 1/2. The Schrodinger equation for (p-^{r±) and 
4'i{'>'i-) therefore becomes. 



-Hm Xr [d/dr^ + (to + 1) /rj_] 

Xr {-d/dr± + m/r±) Hm+i 



E„ 



r 



(3) 



where = -[ti^ / {2M)][d / dr + (l/rj.)d/drx - 

TO^/rj_^] -\- M Lo^r / 2 is the 2D harmonic oscillator. We 
have denoted the energy level as Enm, with n — (0, 1, 2...) 
being the good quantum number in the transverse (ra- 
dial) direction. Each energy level is two-fold degenerate, 
as a result of the time-reversal symmetry satisfied by the 
single-particle model Hamiltonian (Kramer's degener- 
acy). Any state (/>(rj_) = [4'^{Y±_),(j)^{v±_)y is degenerate 
with its time-reversal partner T4'{r±) = {iayC)(j){rj_) = 
[0|(r_L), —(j)'^{r±)]'^, where C is the complex conjugate op- 
erator. Therefore, we may restrict the quantum number 
TO to be non-negative integers, as a state with negative 
TO can always be treated as the time-reversal partner of a 
state with m > 0. To solve numerically the single-particle 
spectrum, we expand the wave-function using the basis 
of 2D harmonic oscillator, 

00 

'^t('^-L) = ^A^kRknrir±), (4) 
00 

(l)i{r±) = ^AikRkm+i{r±) , (5) 



k=0 



where 



Rkr. 



2k\ 



{k+\m\)\ 



1/2 



^1 



2 



is the radial wave-function of "Hm with energy {2k-\- |to 

-fe 



(6) 



l)hio± and /Cj!"' is the associated Legendre polynomial. 



This leads to the following secular equation. 



A, 



— En 



A-^ 
A, 



(7) 



where the vectors A^ and A^ denote collectively the ex- 
panding coefficients {A^fe} and {^4^^}, and the elements 
of matrices (to > 0) are given by, 

(^m)fcfc' = {2k + m + l)Skk'fioj±, (8) 
Mkk' = Xso (Vk' + TO + ISkk' + v^Skk'-ij ^wl9) 

Diagonalization of the secular matrix Eq. ([7]) leads 
to the single-particle spectrum and single-particle wave- 
functions. In numerical calculations, it is necessary to 
truncate the radial quantum number k of the 2D har- 
monic oscillator, by restricting k < fcmax- For Ago ^ 20, 
we find that k^ax — 256 is already sufficiently large 
to have an accurate energy spectrum. With this cut- 
off, the dimension of the secular matrix in Eq. ((T]) is 



2kr. 



512. 



In Fig. la, we present the single-particle energy spec- 
trum at Xso — 5. The ground state single-particle energy 
is plotted in Fig. lb as a function of the dimensionless 
SO coupling constant. In reference to the semi-classical 
zero-point energy E^'^ = —{X'gQ/2)Hui±, the ground state 
energy decrease from huj± to Hu}±/2, when the Rashba 
SO coupling strength Xso becomes sufficiently large. In 
that limit (i.e., Xso > 5), the low- lying energy spectrum 
becomes fairly flat, with a dispersion that is well approx- 



3 




Qi I 0.51 ' ' ■ • ■ ' ■ I 

-20 -10 10 20 1 2 3 4 



so 

Figure 1: (color online), (a) Single- particle energy spectrum 
E„rn = i?nm + (A|o/2)/ki;^ at Xso = 5, measured in reference 
to the semi-classical zero-point energy —{\%Q/2)huj±. (b) 
Ground state single-particle energy, Eoo = -E'oo + (A|o/2)/kJx, 
as a function of the dimensionless Rashba SO coupling con- 
stant. The energy is plotted in units of huj±. 



and 

P2D {E)^2 ^ 5{E-Enmn^), (16) 

in 2D and 3D, respectively. 

For given small numbers of particles N , we can calcu- 
late the low-lying energy levels and then sum explicitly 
the number equations, Eqs. ([12]) and ((T3l) . Once the 
chemical potential is determined at a given temperature 
T, we calculate the occupation of the ground state. 



exp [(£;o - m) AsT] - 1 ' ' ' 

where the single-particle ground state energy Eq — Eqq 
in 2D and Eq = Eqq + hajj2 in 3D. The BEC transi- 
tion temperature Tc can be determined from dPN^/dT'^, 
which exhibits a maximum at Tc (29| . 

III. SEMI-CLASSICAL DENSITY OF STATES 



imated by pl.[l4j. 

-^nm — 



Xso _|_ _|_ !■) ^ m (m + 1) 
2 2 2A|o 



hujA_. (10) 



In 3D, because the motions in xy-plane and z-direction 
are decoupled, the single-particle energy spectrum is 
given by. 



En 



— En 



(11) 



where = 0,1,2... is a good quantum number for the 
axial motion. 

At finite temperature T, the total number of particles 
is given, in the grand-canonical ensembles, by the sum 



in 2D and by the sum 



N 



E 



71,171.71^—0 



exp [{E„ 



fj.) /ksT] - 1 



(12) 



(13) 



in 3D, where /i is the chemical potential and the factor 
of 2 arises from the Kramer's degeneracy. The sum can 
be rewritten as an integral over the energy, in the unified 
form, 



N 



+ 00 



dE- 



p{E) 



exp [{E - ix) /ksT] - 1 ' 
with the DOS p {E) given by 



P2D 



[E) = 2 J2 HE -En 



(14) 



(15) 



For large numbers of particles, it is useful to consider 
a semi-classical approximation by using continuous en- 
ergy spectrum [^]. The level spacing, typical of hujj_ or 
huJzj is assumed to be negligibly small, compared with 
the thermal energy fc^T. Thus, the relevant excitation 
energies, contributing to the sum in Eqs. (IT^ and (IT51) . 
are much larger than the level spacing. The accuracy 
of the semi-classical approximation can be tested a pos- 
teriori by comparing the semi-classical result with the 
numerical discrete summation. 



A. 2D density of states 



In 2D, the semi-classical DOS can be written as, 
,.=±i (27r) 



where E]^siri_) = h'^k\/ {2M) + sXRk^ + Muj\r\/2 is the 
semi-classical energy in phase space (r , ) . Because of 
the Rashba SO coupling, the semi-classical energy splits 
into two helicity branches as indicated by s = ± (see 
Appendix A). By integrating out the spatial degree of 
freedom, we obtain that. 



fy^i-p2D (^) = E / ^-Ldfc^e 



k±_ + s\ 



SO 



E 

huj I 



(19) 

where kj_ = k±a±_ is the dimensionless wave-vector, 
E = E + (A|o/2)fiw_L is the energy measured in ref- 
erence to the semi-classical zero-point energy Eq'^ = 
— {X'^Q/2)huj±, and 8(-) is the Heaviside step function. 
The integration over the wave-vector can be calculated 
explicitly as well. We finally arrive at. 
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r 0, {E< E^^) ; 

f>^±P2D (E) = <^ 2Xso [2E/ (hiu^) + Xlo] ^'^ , (E^- <E<0); (20) 
[2E/{hLo^) + 2Xlo, iE>0). 

In the absence of Rashba SO coupling (Aso = 0), we recover the usual expression for the 2D DOS in harmonic traps, 
P2D (E) = 2E/{hw±)'^<d{E), for a two-component system Q. 



B. 3D density of states 



Likewise, we calculate the semi-classical DOS in 3D, which is given by, 

pTDiE)^Y. I ^5[EUr)~E]^ 



(21) 



where the semi-classical energy now takes the form i?ks(r) = h^k\/ {2M) -\- s\Rk:i_ + h^k'l/ {2M) + M{ijj\r']_+uj'lz'^) /2. 
The integration over r and kz can be done by introducing a new variable = fi^kl/ {2M) + M{uj\r\ + a;^2;^)/2 and 
by converting the variables of integration drdk to dtdkj^. This leads to. 



oo 

hw^pll {E) = Y, [ k±dk 



E 



k± + sX 



so 



6 



E 



k± + sX 



so 



(22) 



By explicitly integrating out k±, we obtain. 



0, (E < E^-) ; 

h^zP^D (E) = { (2Aso/3) [2E/ (fu^^) + Xl^f^ , (E^- < E < 0); 

[i?/(fi^^) + A|o]'-A|oA {E>0). 



(23) 



In the absence of Rashba SO coupling, we recover the ex- 
pression pfjj {E) — E^/ (^h'^Lj^uJz) &{E) for 3D harmonic 
traps 

It is easy to check that the 2D and 3D DOS is related 
by 



dpTo {E) 
dE 



Pfo {E) . 



(24) 



This is due to the decoupled motion in the x — y plane 
and z direction, which leads to the observation that the 
3D energy spectrum may alternatively be viewed as a 
collection of 2D spectra with regular spacing hujz- 



use 



S{x) = 2fs{x; r = hoj^) - fs{x; T = 2nu^), (25) 

as an extrapolation to the zero-broadening limit (F = 
0). We find that the semi-classical expressions for DOS, 
Eqs. (fT9|) and ([22]), works extremely well over a very 
broad range for energy. The most significant discrepancy 
occurs at the lowest energy level, E ~ —{X'gQ/2)hu!±, as 
anticipated. 



IV. CRITICAL TEMPERATURE AND 
CONDENSATE FRACTION 



C. Test of the semi-classical DOS 

In Fig. 2, we compare the semi-classical 2D and 3D 
DOS with these obtained by summing over the discrete 
single-particle energy spectrum using Eqs. (jl5p and 
(|16p . In the numerical summation, we simulate the 6- 
function 6{x) by a Lorentzen line shape with broadening 
F, f5{x;T) = (F/7r)/(x2 -|- F^). Roughly, the resulting 
DOS depends linearly on F at F ^ hujj_. Therefore, we 



We are now ready to calculate the critical tempera- 
ture and condensate fraction for large number of parti- 
cles. With the semi-classical DOS p""^ (E), the number of 
particles could be rewritten as [4], 

^-^0+ dE / \ / (26) 

Ji5- exp [{E - p) /ksT] - 1 

where the ground state population iVo is singled out and 
the finite sum over the excited states in Eqs. (jl2p and 
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Figure 2: (color online). The semi-classical density of states 
in 2D (a) and 3D (b), in units of l/{f\hj±) and l/(hwz) re- 
spectively, are shown as a function oi E = E + (A|o/2)^^ 
at different Rashba SO couplings (solid lines). The energy 
E is in units of tujj^. For comparison, the symbols plot the 
results obtained by the numerical summation, see, Eqs. (|15p 
and p6|l . The simulation of the delta- function is described in 
the text. 



(|13p is replaced by an integral. Accordingly, we have set 
the lower-bound of the integral to be the semi-classical 
zero-point energy Eq'^ = —{\'gQ/2)huj±. When BEC oc- 
curSjthe chemical potential approaches to Eq'^ from be- 
low [J| . The critical temperature Tc is determined by the 
condition, 



N 



+ 00 



dE- 







exp 


E/keT, 


- 1 



(27) 



where E = E — Eq'^, and the condensate fraction at T < 
Tc can be calculated by. 



No 1 ~ P''(e + E^ 

— = 1 / dE P : 

^ ^ Jo exp E/ksT 



- 1 



(28) 



As we shall see, these equations can be conveniently 
solved by introducing e = E / {ksT) and 



a{T) = \so 



(29) 



2D 



In 2D, the equations for the critical temperature and 
condensate fraction becomes. 



N = 



and 



No 
N 



= 1 - 



I2D [a (Tc)] 



I2D [a (T)] 
I2D [a (Tc)] ^ 



(30) 



(31) 



respectively. Here the integral I2D [a] takes the form, 

i,„H^r^.^¥^, (32) 



where the dimensionless DOS pj^j (e; ct) is given by. 



^20 (e; a) 



0, 

2aV2e, 
2e + a^ 



(6 < 0) ; 

(0 < e < qV2); 

(6>«V2). 



(33) 



Therefore, 



I2D [a] = V^TTaC ( - 



+ 00 



de- 



(34) 



Here C(-) is the Riemann C function. X2D [a] depends 
implicitly on the temperature through the dimensionless 
parameter a(r). It is clear from Eq. (|29|1 that for a 
given SO coupling A50, the dimensionless parameter a 
at the critical temperature Tc always scales to zero in the 
thermodynamic limit iV — >■ 00. This is understandable as 
a finite SO interaction modifies only the low-lying energy 
states, whose occupation becomes negligible as iV — )■ cx). 

In the absence of SO coupling, X2d[q; = 0] = 2(^(2) = 
7r^/3, we recover the standard results in 2D J^], 



fes^io) (\so = 0) = i (37V)^/' 



(35) 



and No/N = 1 — {T /Tc'^^Y . Here, we use the superscript 
"0" to indicate the semi-classical result. For a large SO 
coupling, I2D [a ^ 1] = ^/2TTaC, (3/2), we find 



fesTjo) (Aso » 1) = 



1 



(27r) 



1/3 



N 



AsoC(3/2) 



2/3 



(36) 

and No/N = 1 - [T/T^^^ff^. Thus, for a given number 
of particles, with increasing SO coupling the dependence 
of 2D critical temperature on the number of particles 
changes from N^^'^ to N'^/'^. Using a ^ 1, the strong- 
coupling limit is reached when 



Aso » (2^)-'/' 



N 
C(3/2) 



1/4 



(37) 
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Figure 3: (color online). The integrals T2D and I315 as a func- 
tion of the dimensionless parameter a = \so[h^ -L / {ksT)]^^'^ 
(symbols) . The solid lines show the empirical fit which agrees 
numerically within 0.5% relative error (see the text for the 
empirical formalism). 
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Figure 4: (color online). 2D critical temperature as a func- 
tion of the SO coupling at different numbers of particles, as 
indicated. The dot-dashed lines show the limiting behavior 
at large SO coupling, Eq. (l36)l . 




Figure 5: (color online). 2D condensate fraction at Aso = 5 
and at different numbers of particles. For the case of A'^ = 10''', 
the dashed and dot-dashed lines show respectively the strong- 
coupling and zero-coupling result, Nq/N — 1 — {T/Tc^^)"^ and 
No/N ^1-{T/T^''''f'^. 



In Fig. 3, we show I2D [a] as a function of the dimen- 
sionless parameter a. Empirically, we find that I2D [c«] — 
\/2^aC (3/2) -I- 2C(2)e-i-84"-o i3"'^ within 0.5% relative 
error. Fig. 4 reports the critical temperature as a func- 
tion of SO coupling at several numbers of particles (solid 
lines). It decreases significantly at moderate SO coupling 
(A50 ~ 10) and number of particles (i.e., ^ 10*). The 
strong-coupling results Eq. ([55)) have also been plotted 
using dot-dashed lines. Finally, in Fig. 5, we present the 
condensate fraction at A50 — 5 and N — lO'^, 10^, and 

CX). 



B. 3D 



In 3D, similarly we obtain that 



and 



TV 



^_.T\'l,u [a{T)] 



(38) 



(39) 



where A = u}z/u}± is the aspect ratio of the harmonic 
trap, the integral I^u [a] is given by, 



^3D \a\ = 



e^-l ' 

and the dimensionless DOS p^'jy (e; a) is, 

r 0, (e < 0) ; 

pTd (e; «) = { (4y2a/3)e3/2, (0 < e < aV2); 

I e2 + ea2-aVl2, {e>a^/2). 



(40) 



(41) 



Explicitly, we find that 

+00 

T3D H = v^aC (0 + / d^^^^ (42) 



where h{e) = + ea^ - - {AV2a/3)e^/^ . We 

plot [a] in Fig. 3, together with an empirical fit, 
I3D [a] = V2^aC (5/2) + 2C(3)e-i-40"-"-3""'. 

At Aso — where lania = 0] = 2^(3), we obtain 



AsTjo) (Aso = 0) 



A^A 



2C(3) 



1/3 



(43) 



and No/N = 1-{T/tP)^, recovering the well-known 3D 
result for a trapped spin- 1/2 Bose gas In the limit of 
large SO coupling where I3D [a » 1] = V^na^ (5/2), we 
find instead 

2/5 



fesTW (Aso » 1; 



{2n 



,1/5 



iVA 
AsoC(5/2) 



(44) 
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and No/N = 1 - (T/T^^'^f''^. Thus, for given iV, with 
increasing SO couphng the power-law dependence of 3D 
critical temperature on number of particles changes from 
7V^/^ to iV^/^. We estimate that the strong-coupling re- 
sult is applicable if 



\so » (2^)-^/^' 



1/6 



(45) 
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Figure 6: (color online). 3D critical temperature as a function 
of the SO coupling strength at A'' = 10^ and A*' = 1.8 x 
10"". The aspect ratio of the harmonic trap is A = cj^/ojx = 
\/8. The dot-dashed lines show the critical temperature in 
the strong-coupling limit, Eq. (144(1 . The inset reports the 
condensate fraction at \so ~ 5. 

In Fig. 6, we report the effect of the SO coupling 
on 3D critical temperature. To make a connection with 
the NIST experiment we have used a realistic aspect 
ratio of the trapping potential and number of particles, 
A = VS and iV = 1.8 x 10^. We also consider the case 
with a small number of particles = 10^. At the typical 
SO coupling strength Ago 10 [H, the reduction of the 
critical temperature is about 20%, which is in reach of 
current experiments. The inset shows the condensate 
fraction at Xso — 5. 



V. FINITE SIZE CORRECTION TO Tc IN 3D 

We now turn to consider the finite size correction to the 
semi-classical results, which arises from the discreteness 
of the single-particle energy spectrum ^26i . j27i] . The semi- 
classical results are obtained using the semi-classical ap- 
proximation for the excited states and setting the chemi- 
cal potential to the semi-classical zero-point energy Eq'^. 
To the leading order, the finite size correction can be in- 
cluded by still employing the semi-classical description 
for the excited states, while keeping the quantum value 
H — Eo for the chemical potential at the transition (2^ . 
Here, Eq > Eq'^ is the single-particle energy of the ground 



state. It is Eqq in 2D and Eqq + huJz/2 in 3D; see, for ex- 
ample. Fig. (lb) for Eqq as a function of the SO coupling 
strength. The discreteness of the excited energy spec- 
trum gives rise to higher-order finite size corrections. In 
the following, we focus on the finite size correction to the 
3D critical temperature. 

Using the quantum value fj, = Eq for the chemical po- 
tential, the 3D critical temperature is determined by. 



A^ 



+ 00 



dE- 



pTd (E) 



E, exp [{E - £^0) IkeTc] - 1 ' 

+00 (e + E^- 

dE 

exp E/ksTc - 1 



AEj 



(46) 
(47) 







where in the second line we have introduced E — E ^ E[ 
and AE = Eo- E^" > 0. Compared with Eq. (I?7l) . 
the 3D DOS is slightly up-shifted by an amount AE. As 
AE ~ huj± is the smallest energy scale, using Eq. ([24|) 
we may write pl%{E + E^" + AE) ~ p's%{E + E^") + 
{AE/f}Wz)p^E){E + i?o^). Therefore, using the integrals 
T2D and the equation for the critical temperature is 
given by 



A^A 



huj I 



3 r 



'-ZD 



[«m)] 



AE 



I2D [a {Tc 



(48) 

In the absence of the SO coupling, I2D — 2C(2), I3D — 
2C(3) and AE = hoj± + hujz/'^, it is easy to verify that 
the transition temperature Tc is given by the law. 



^^1 

TP " 



C(2) 



[2C(3) 



2/3 



{2huj± + huj^) /3 



(w^Wz) 



1/3 



(49) 



which is known in the literature |25l - l27| . 

In Fig. 7, we report the 3D transition temperature 
with the leading finite size correction, as shown by dashed 
lines. We find a sizable correction at small number of 
particles (i.e., A^ = 10"^). For the experimentally realistic 
number of particles, i.e., A^ = 1.8 x 10^, however, the cor- 
rection becomes mild. As a benchmark to our analytic 
treatment for Tc, we also show by symbols the critical 
temperature for small number of particles, calculated by 
the discrete sum for the ground state population Nq, Eq. 
pri) . At relatively small SO coupling (i.e.. Ago < 5), our 
analytic treatment works very well. However, for large 
SO coupling, the single-particle level splitting between 
the ground state and the first excited state becomes in- 
creasing small. We then may have to take into account 
the discreteness of the low-lying excited energy levels. 



VI. CONCLUSIONS 

In summary, we have investigated the critical temper- 
ature and condensate fraction of a harmonically trapped 
ideal Bose gas in the presence of Rashba spin-orbit cou- 
pling, by using either the exact numerical summation for 
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Figure 7: (color online). 3D transition temperature as a 
function of the SO coupling. Tiie solid lines show the semi- 
classical predictions and the dashed line gives the results with 
inclusion of the leading finite size correction. The empty cir- 
cles are calculated using the numerical summation for A'o with 
discrete energy spectrum, i.e. Eq. (|17ll . The critical temper- 
ature is then determined from the peak position of (fNo/dT^ 



small number of particles or the analytic semi-classical 
approach for large number of particles. The leading fi- 
nite size correction to the semi-classical approximation 
has also been considered. We have found pronounced ef- 
fect of the Rashba SO coupling. For the experimentally 
realistic number of particles [N ^ 10^) [1:], the critical 
temperature is reduced by more than 20% in magnitude 
at a moderate SO coupling. This reduction is readily ob- 



servable in current experiments. Moreover, in the limit 
of strong SO coupling, the critical temperature scales as 
N^^^ and A^^/-^ in three and two dimensions, respectively, 
which should be contrasted with the scaling law of N^^^ 
and N^/^ in the absence SO coupling. Our investigation 
of critical temperature can be easily extended to include 
a weak repulsive interaction, by using mean-field Hartree- 
Fock theory 
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Appendix A: Density of states of a 3D homogeneous 
SO coupled system 

In free space, the single-particle Hamiltonian with 
Rashba SO coupling. 



has the dispersion, 



-iXnidy + idx) 



Evs — 



2M 



2M 



(Al) 



(A2) 



Here s — ± denotes the two helicity branches. The DOS, 
given by piE) = (l/V) Ek[SiEk+ ~ E) + 5{E^^ - E% 
can be calculated analytically. We find that. 



p{E) 



0, 

7r/2, 



yj2E/EB, + 7r/2 - arctan yj2E/EB, 



{E < -Er/2) ; 
{-Er/2 <E<0); 

iE>0). 



(A3) 



where Er = AlXj^/h"^ is the characteristic energy related to the SO coupling. This result was reported by Hui Zhai 
in Ref. [6| (see for example, their Fig. 2b). By introducing a Fermi wave-vector kp = (Btt^A^/F)^/^, Fermi energy 
Ep = fi^kp/ {2AI), and dimensionless SO coupling strength Ae// = APXR/{h'^kp), the DOS can be written as. 



P{E) 




arctan 



^E/{xl,,Ep) 



^2 <E<Q); 
{E>0). 



(A4) 



We show in Fig. 8 the DOS at different SO coupling 
strengths. 
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